Tomographic reconstruction of molecular orbitals with twofold mirror antisymmetry: 

overcoming the nodal plane problem 
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We propose a new method to overcome the nodal plane problem for the tomographic reconstruc- 
tion of molecular orbitals with twofold mirror antisymmetry in the length form based on high-order 
harmonic generation. It is shown that, by carrying out the reconstruction procedure in the rotat- 
ing laboratory frame using the component of the dipole moment parallel to the electron recoUision 
direction, the nodal plane problem is avoided and the target orbital is successfully reconstructed. 
Moreover, it is found that, the proposed method can completely avoid the additional artificial lobes 
found in the results from the traditional method in the velocity form and therefore provides a more 
reliable reproduction of the target orbital. 
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PACS numbers: 32.80.Rm, 42.65.Ky 

The fast progress in strong-field physics has made it pos- 
sible to observe the structure and ultrafast electron dy- 
namics in atoms and molecules with Angstrom and at- 
tosecond resolutions [THS]. In the recent decade, a fas- 
cinating application known as the molecular orbital to- 
mography (MOT) based on high-order harmonic gener- 
ation (HHG) is raised (THIS]. With the MOT scheme, 
a two-dimensional projection of the molecular orbital on 
the plane orthogonal to the pulse propagation direction 
can be imaged after measuring the HHG spectra from 
different alignments of the molecule. In detail, by cali- 
brating the measured HHG signal of the target molecule 
with that of a reference system with the same ionization 
energy, the recombination dipole matrix element of the 
molecular orbital is extracted. Then based on the Fourier 
slice theorem, the target molecular orbital can be recon- 
structed by performing the inverse Fourier transform of 
the recombination dipole matrix elements (see Refs. [M]- 
[T6] for more information). 

The reconstruction can be in principle performed in 
both the length form (LF) and the velocity form (VF). 
However, the applicabilities of the two forms are differ- 
ent. It is generally suggested to apply the velocity form if 
the target orbital has one or more nodal planes contain- 
ing the coordinate axes in the molecular frame, because 
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one would run into numerical problems in the length form 
p!4HT6] . In this work, we show that this problem can be 
avoided if the orbital is antisymmetric with one nodal 
plane, and the problem only holds for the orbitals with 
twofold mirror antisymmetry. Accordingly, we propose 
a new method to overcome the nodal plane problem in 
reconstructing this kind of orbitals in the length form. 
Taking a iVg orbital as an example, we compare the recon- 
structed result with that obtained from the traditional 
method in the velocity form, which shows that the pro- 
posed method provides a more reliable reproduction of 
the target orbital. Finally, we discuss the scope of appli- 
cation of this method and show that it can be extended 
to more types of twofold antisymmetric orbitals to over- 
come the nodal plane problem. 

The scheme of MOT is illustrated in Fig. 1. (xyz) is 
the molecular frame and (x'y'z') is the laboratory frame. 
The probe laser pulse propagates along the z' axis and 
is linearly polarized parallel to the x' axis, along which 
the continuum electron recollides with momentum k to 
generate the HHG radiation. is the alignment angle. 

In the length form, the traditional reconstruction 
method can be expressed as p!4Hl6] : 

^^^{x,y) 

where Tk^r denotes the inverse Fourier transform. Since 
the reconstruction is derived in the molecular frame, the 
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FIG. 1: (Color online) Illustration of the MOT scheme: (xyz) 
is the molecular frame and (x'y'z') is the laboratory frame. 
The probe laser pulse propagates along the z axis and is 
linearly polarized parallel to the x axis, along which the con- 
tinuum electron recollides with momentum k to generate the 
HHG radiation. Q is the alignment angle of the molecule. 

X and y components of the dipole moment (d^ and dy) 
are required. In experiments, the HHG radiation in both 
the x' and y' directions are measured first to extract 
G^ll and d\_ using a polarizer formed by a pair of silver 
mirrors [HI [17] , then d^ and dy are calculated by pro- 
jecting (i|| and d^ onto the x and y axes respectively: 
dx = d\\ cos{6) — d±sm{6), dy = d\\ sm{6) + dj_cos{6). 
In the Fourier space, kx = /ccos(6>), ky = ksm{0) and 
k'^/2 = uj according to the energy conservation, where 
uj is the harmonic frequency. Although the same or- 
bital can in principle be reconstructed from both the 
X and y components of the dipole moments, they will 
most likely not give the same result due to the limited 
discrete sampling in Fourier space. Therefore, the final 
reconstructed orbital is generally defined as ^^y{x^y) = 
^[^^^{x^y) + ^y^{x^y)] to hopefully average out the 
distortions. 

In Eqs. (1) and (2), if the target orbital has a nodal 
plane containing the y or x axis, one will run into numeri- 
cal problems when dividing hy x or y respectively [T4HT6] . 
In detail, if the orbital has a nodal plane containing the 
y axis for example (i.e. is antisymmetric with respect to 
the y axis), the result of the inverse Fourier transform in 
the sampled spectral range in Eq. (1) Tk^r[dx] equals 
a nonzero value at x=0 and one will confront the nu- 
merical problem when performing the division by x since 
^^j^ k^r[ x\ _ ^^ Q^ ^^^ other hand, if the orbital is 

symmetric with respect to the y axis (i.e. without the 
nodal plane), Tk^r[dx] -^ and the division ^^"^^ ""^ 
tends to a finite value for x -^ 0, so one will not run 



into the nodal plane problem. These results are due to 
the symmetry property of the orbitals and the feature 
of the inverse Fourier transform in the limited k-space. 
Besides the above discussed orbitals with symmetry or 
antisymmetry, the orbital can also be nonsymmetrically 
distributed, the reconstruction of which was discussed in 

m- 

In previous works, it is generally suggested that one 
should employ the velocity form to avoid the nodal plane 
problem. Actually, this problem can also be avoided by 
choosing appropriate component of the dipole moment in 
the length form, if the orbital has only one nodal plane. 
Take the highest occupied molecular orbital (HOMO) of 
acetylene for example [19] , which is a tTu orbital with only 
one nodal plane containing the x axis as shown in Fig. 
2(a). The target orbital can be reconstructed by using 
only dx (according to Eq. (1)) as shown in Fig. 2(b), 
while can not be reconstructed using dy (according to 
Eq. (2)) as shown in Fig. 2(c). The unexpected structure 
observed along the x axis in Fig. 2(c) is just owing to the 
nodal plane problem. In Fig. 2(d), we also present the 
reconstructed result by using the traditional method in 
the velocity form [HJ fl4| . Comparing panel (b) with (d) , 
it is found that the target orbital is better reconstructed 
in the length form. In the simulation, the dipole moment 
is calculated by d = (^|r|e*^'^), where ^ is the Hartree- 
Fock (HF) orbital obtained from the Gaussian 03 code 
[2Q1 . Throughout this work, the spectrum is sampled 
by the odd harmonics of the 800-nm laser field in the 
range of harmonics 17-51 [7 . The alignment angle is 
scanned from 0° to 350° in step of 10°. We assume a 
perfect alignment of the molecules in the simulation and 
the molecules are considered to be fixed in space during 
the interaction. 

However, when the target orbital has two nodal planes 
containing the x and y axes respectively {i.e. with the 
twofold mirror antisymmetry), one will run into the nu- 
merical problem in both Eqs. (1)(2) and the nodal plane 
problem remains unavoidable. Take the HOMO of CO2 
with TTg symmetry for example [11 . Fig. 3(a) shows the 
HF orbital obtained from the Gaussian 03 code, in which 
two nodal planes divide the orbital into four lobes with al- 
ternating signs. Fig. 3(b) shows the reconstructed result 
with the traditional method in the length form. What 
one can see is only the cross-shaped structure resulting 
from the nodal plane problem. 

To overcome this problem, we propose a new recon- 
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FIG. 2: (Color online) 2D projections of the HOMO of acety- 
lene obtained (a) from HF calculation, (b) according to Eq. 
(1), (c) according to Eq. (2), and (d) with the traditional 
method in VF. 
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FIG. 3: (Color onhne) 2D projections of the HOMO of CO2 
obtained (a) from HF calculation, (b) with the traditional 
method in LF, (c) with the proposed method, and (d) with 
the traditional method in VF. 



structing method in the length form. The basic idea is 
to reconstruct the orbital using d\\ directly in the labo- 
ratory frame instead of projecting d\\ and d_\_ at different 
onto the fixed x and y axes in the molecular frame. 
According to d\\{k,0) = {'^\r\\^Q\e'^^^'^), the target orbital 
is reconstructed following 



*f(x,y) = ^ 



J^k- 



9{ka:,ky)] 



^\\,e{x,y) 



(3) 



For any point P(x,y), r\\^0{x^y) = x cos{6) -\-y sm{6) is the 
projection of the vector OP onto the x' axis (as shown 
in Fig. 1). In this method, the division is performed in 
the rotating direction ry 51 (signal along ry 51 divided by 
the coordinate on ry 51) rather than in the fixed x and 
y directions perpendicular to the nodal planes. In this 



frame, the inverse Fourier transform can be expressed as 

Tk^r[d\\^e] = d\\^0 ex.-p[ik{x cos -\- y sin 0)]dk (4) 

= d\\^QCOs{kr\\^Q)dk -\- i I d\\^Qsiii{kr\\^Q)dk (5) 

To explain why the proposed method can avoid the 
nodal plane problem in a concise form, we assume that 
the molecular orbital is real- valued. (An arbitrary com- 
plex phase of the orbital does not change the result of the 
reconstruction and can actually be set to zero by apply- 
ing a global phase shift 0|T5].) In this case, if the target 
orbital ^ is symmetric about the origin of the molecular 
frame, d\\^Q is purely imaginary- valued [9 . Therefore, the 
reconstructed orbital should be 
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According to Eq. (7), since the limitation of the divi- 



sion lim 



id||,g,fcsin(fer||,g) _ 



ikd\\ Q f^ equals a finite value 



for each alignment angle ^, the nodal plane problem is 
avoided. 

The reconstructed result with this proposed method 
for CO2 is shown in Fig. 3(c). The main structure of 
the HOMO of CO2 is successfully reproduced, and the 
main error is only the diffusion of the orbital in the outer 
region. As a comparison, we also present the result recon- 
structed by the traditional method in the velocity form in 
Fig. 3(d). In this figure, one could see eight or more lobes 
in the reconstructed result (the same as in [rP, T!4]). The 
additional artificial lobes would be misleading for peo- 
ple to learn the structure of an "unknown" orbital. This 
problem is also found in Fig. 2(d): the reconstructed or- 
bital from the traditional method in the velocity form is 
seriously distorted by the additional lobes. 

To compare the results more clearly, the slices of the 
HF orbital, the reconstructed orbitals ^m"^ and ^^^ of 
the HOMO of CO2 are presented in Fig. 4. We compare 
two groups of slices. Panels (a) and (b) show the cuts 
along the two diagonal directions respectively, and panels 
(c) and (d) depict the cuts along y = 1 a.u. and x = 2 a.u. 
respectively. The comparison shows that, although the 
slices of ^^^ match better with those of the HF orbital 
for the amplitudes of the peaks, they have serious oscil- 
lations which do not exist in the slices of the HF orbital. 



These oscillations just correspond to the additional lobes. 
On the other hand, regarding ^m'^ reconstructed by the 
proposed method, although the result does not match the 
HF orbital exactly for the absolute values and converge 
to zero more slowly, it completely avoids the oscillations 
{i.e. the additional lobes). The comparison indicates 
that, without the misleading additional lobes, the pro- 
posed method provides a more reliable reproduction of 
the target orbital. It is also worth noting that, the qual- 
ities of the reconstructed results are also affected by the 
choice of the sampled spectral range for both the tradi- 
tional method in velocity form and the proposed method. 
Since the width of the detected spectrum is limited in ex- 
periments, the sampled spectral range should be carefully 
chosen. This has been systematically discussed in [M], 
and it is summarized that the spectral range should con- 
tain a characteristic spatial frequency of the orbital and 
a "well balanced" amount of the positive and negative 
spectra amplitudes. In this work, the reconstructed re- 
sults are expected to be better for both reconstructing 
methods, if the spectral range is further extended. How- 
ever, one might confront more experimental problems to 
achieve the broader detected spectral range. 
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FIG. 4: (Color online) Comparison of the slices of the HF 
orbital, the reconstructed orbitals ^|^ and ^^^ for CO2: 
(a)(b) along the diagonal directions, and (c)(d) along two 
straight lines parallel to the x and y axes respectively. 



Before conclusion, we would like to briefly discuss the 
scope of application of the proposed method. Although 
the applicability of this method is not limited by the 
nodal planes of the target orbital, it can not be used for 
all kinds of orbitals either. It requires the target orbital 
to be symmetric about the origin of the molecular frame, 
or one will again run into numerical problems when divid- 
ing by ry according to Eqs. (5-7) and related discussion. 
For the orbitals with twofold antisymmetry, which have 
the nodal plane problem, this requirement is satisfied. 
Therefore, this method can be applied to various types 
of molecular orbitals with the twofold antisymmetry to 
overcome the nodal plane problem. 



In summary, the nodal plane problem can be avoided 
for the antisymmetric orbitals with one nodal plane, and 
only holds for the orbitals with twofold antisymmetry for 
the traditional method in the length form. To overcome 
this problem, we propose a new method by performing 
the reconstruction procedure in the rotating laboratory 
frame using the component d\\ parallel to the electron 
recollision direction. With this method, the HOMO of 
CO2 is successfully reconstructed and we further com- 
pare the reconstructed result with that from the tra- 
ditional method in the velocity form. The comparison 
shows that, our proposed method completely avoids the 
additional artificial lobes found in the orbitals obtained 
in the velocity form and therefore provides a more reli- 
able reproduction of the target orbital. Finally, it is dis- 
cussed that this method can be extended to more types 
of orbitals with the twofold antisymmetry. 
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